Dynamical Decoupling of Spin-Clusters using Solid State NMR 
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In this work we experimentally study the efficiency of various dynamical decoupling sequences for 
suppressing decoherence of single as well as multiple quantum coherences on large spin-clusters. The 
system involves crystallites of a powdered sample containing a large number of molecular protons 
interacting via long-range dipole-dipole interaction. The multiple quantum coherences are prepared 
by progressively creating correlations in the spin lattice using standard pulse sequences implement- 
ing two-quantum average Hamiltonian. The spin system is then subjected to various dynamical 
decoupling sequences, followed by conversion into observable single quantum coherence by using 
time-reversal sequence. The experiments reveal superior performance of the recently introduced 
RUDD sequences in suppressing the decoherence. 
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I. INTRODUCTION 

The study of dynamics and control of quantum many 
body systems has renewed interest in the field of quan- 
tum information. While encoding information onto a 
quantum channel can potentially speed up certain com- 
putations and allow secure data transmission, the prac- 
tical realization of these applications are hindered by the 
extreme sensitivity of the quantum channel to environ- 
mental noises. Systems based on nuclear spin-clusters is 
one among the various architectures being investigated to 
realize quantum channels. Several experimental demon- 
strations of quantum information processing (QIP) using 
solid-state nuclear magnetic resonance (SSNMR) have 
already been reported [IHll- By sophisticated control 
of spin-dynamics it is in principle possible to achieve a 
larger number of quantum bits (qubits) using SSNMR, 
because of the availability of large spin-clusters coupled 
mutually through long-range dipole-dipole interactions. 
However in such a spin-cluster, fiuctuating local fields at 
the site of each spin induced by its environment leads to 
the decoherence of the encoded quantum information. 

Due to the availability of large spin-clusters it is possi- 
ble to prepare coherences of large quantum numbers by a 
widening network of correlated spins evolving under two- 
quantum average Hamiltonian [g, 0] ■ These higher order 
coherences are not directly observable as macroscopic 
magnetizations, but can be converted into observable 
single quantum coherence (SQC) using a time-reversed 
two-quantum average Hamiltonian. This method, often 
known as a 'spin-counting experiment' has been used to 
study the evolution of coherences of large quantum num- 
bers exceeding 4000 . 

Under the standard Zeeman Hamiltonian any spin- 
coherence is a non-equilibrium state and decays via var- 
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ious relaxation processes, ultimately leading to the equi- 
librium longitudinal magnetization. It has long been dis- 
covered that the decay process of the spin coherence can 
be prolonged by applying a series of spin flips at regu- 
lar intervals of time. This sequence known as 'CPMG 
sequence', not only refocuses the effect of spectrometer 
inhomogeneities, but also reverses the phase evolution of 
the coherences under the random fluctuations, provided 
the spin flips are applied sufficiently frequent [ll|, ll^ . 
Such a dynamical method for the suppression of deco- 
herence of a qubit due to its interaction with environ- 
ment is often termed as 'dynamical decoupling' (DD) 
[13} . Recently Uhrig introduced a non-periodic spin-flip 
sequence which he proved theoretically to provide opti- 
mal decoupling performance for dephasing spin-bath in- 
teractions J^J. CPMG and other similar periodic spin- 
flip sequences suppress spin-environment interaction to 
nth order using up to 0(2") pulses, while Uhrig dynam- 
ical decoupling (UDD) suppresses the same using only n 
pulses. Filter function analysis of dynamical decoupling 
indicates that, in a high frequency dominated bath with 
a sharp cutoff, UDD works well provided the frequency 
of the spin- flips exceeds the cutoff frequency fla-KLTf . On 
the other hand when the spectral density of the bath has 
a soft cutoff (such as a broad Gaussian or Lorentzian), 
the CPMG sequence was found to outperform the UDD 
sequence (1J- 25J. Later on, UDD has been generalized to 
suppress simultaneously both transverse dephasing and 
longitudinal relaxation of a qubit [29|. Also, the original 
sequence for UDD is based the assumption of instanta- 
neous spin-flips, which requires inflnite bandwidth. More 
recently, Uhrig provided an improved sequence - 'realistic 
UDD' (RUDD) for practical implementations with finite 
bandwidth [27| . 

Most of the theory and experiments of DD sequences 
are for single spin systems. Du et al have studied DD 
of electron spin coherence in solids [ll] , while Suter and 
co-workers have reported systematic experimental com- 
parisons of various DD schemes on an ensemble of single 
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spins in SSNMR [19|,l25|. 

Later Agarwal has shown using theoretical and nu- 
merical calculations that even entangled states of two- 
qubit systems can be stored more efficiently using UDD 
[28,]. Similar theoretical studies were also carried out by 
Mukhtar at. al. [l^. Experimentally, Wang et. al. have 
studied DD on electron-nuclear spin pairs in a solid-state 
system 30], and Soumya et. al. have studied the perfor- 
mance of UDD an a two-qubit liquid-state NMR system 
p. 

In this article, we report the experimental study of 
performance of various DD schemes on an extended net- 
work of spin- 1/2 nuclei forming a large spin cluster. This 
article is organized as follows. Next section briefly de- 
scribes the method of preparing and detecting multiple 
quantum coherences (MQC) in SSNMR, section III sum- 
marizes the construction of various DD sequences, and 
the experimental details are described in section IV. Fi- 
nally we conclude in section V. 



II. MULTIPLE QUANTUM SSNMR 
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The SSNMR Hamiltonian for a spin cluster with M 
spin- 1/2 nuclei is 
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where the Zeeman and the secular part of dipolar inter- 
action are respectively. 
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|33 | . Here I* and are spin angular momentum operator 
and its z-component corresponding to the i*'' spin, and Wi 
and Dij are the chemical shift and the dipolar coupling 
constants. The equilibrium density matrix for the above 
Hamiltonian corresponds to the longitudinal magnetiza- 
tion expressed as X^i-^l- The density matrices for the 
longitudinal spin order can be expressed using product of 
longitudinal spin operators, eg. /^/^ • • • . The coherences 
are described by the product of transverse (or of trans- 
verse and longitudinal) spin operators, eg. I^I^Ix ' ' ' ■ 
The transverse spin operators can also be expressed in 
terms of raising and lowering operators: Ix — (/++/-)/2 
and ly = — — /_)/2. The difference between the total 
number of raising and lowering operators gives the quan- 
tum number n of a particular coherence. For example, 
operators I-^I^, and describe zero, single, and 

two-quantum coherences respectively [ssj . 

The pulse sequence for preparing and detecting MQC 
is shown in Fig. [Ija-b). The sequence in Fig. [1^ involves 
preparation of MQC, application of DD schemes, free- 
evolution (tl), converting MQC into longitudinal spin or- 
der (mixing), destroying the residual coherences by trans- 
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FIG. 1: The experimental scheme (a) for studying the perfor- 
mance of DD on large spin-clusters and the 8-pulse sequence 
(b) implementing "Hi (a). In (b). A' — 2 A -I- where 
is the duration of each n/2 pulse. The DD schemes are de- 
scribed in (c-e). The phase 4> is set to x for CPMG, UDD, 
and RUDD schemes, while it is alternated between x and —x 
for CPMGp, UDDp, and RUDDp. 



verse relaxation (tj?), followed by detection after convert- 
ing the longitudinal spin order into SQC. The 8-pulse 
sequence in Fig. [T|d corresponds to the two-quantum av- 
erage Hamiltonian 



(3) 



for a = 0. Preparation and mixing parts involve m- 
cycles of the 8-pulse sequence %m{a) and ■Hm(7r/2) 01 . 
Under the preparation sequence, each uncorrelated spin 
gets correlated with its neighboring spins, creating two- 
quantum coherence. Each pair of correlated spins evolve 
under their neighbors, and this growing network of cor- 
related spins in the form of higher order coherences leads 
to large correlated clusters. The possible quantum num- 
bers and the corresponding cluster size increases with the 
number of cycles. Only even quantum coherences are pre- 
pared as shown in Fig. [5J However, MQCs themselves 
have no detectable macroscopic observables. MQCs can 
be converted back to longitudinal spin order by a time- 
reversal sequence 'Hm (7r/2) (mixing), and subsequently 
be converted into detectable SQCs by a Ti/ly pulse. To 
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FIG. 2: Possible quantum numbers (n) for M-spin cluster 
after excitation with several cycles of 8-pulse sequence (Fig. 
[TJj) implementing two-quantum average Hamiltonian (Hm)- 



C. RUDD 

In RUDD and RUDDp, in addition to the delays be- 
tween the pulses, the pulse durations and amplitudes also 
vary, but the overall sequence remains symmetric [27|. 
The pulse durations are given by 



T 



N+1 



(5) 



where T is total duration of the sequence and N is the 
number of pulses. Here Op is a constant and can be de- 
termined by the allowed bandwidth. We choose = , 
and calculated 9p based on the minimum allowed pulse 
duration; 



separate the MQCs, the relative phase a between the 
preparation and mixing is incremented in proportion to 
the evolution time ti (Fig. [T^). Spurious transverse co- 
herences are suppressed by an extended delay tf/. A final 
purge pulse px is used to keep only the x-component. 



III. DD SCHEMES 



A. CPMG and CPMGp 
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The amplitude aj of j pulse is calibrated such that 
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n. Time instants of the center of each pulse 
is same as in equation (|4]). Using these time instants, 
the delays between the pulses can be calculated as ti — 
tn+1 ^h- and r, = tj - tj_i - - t^~^ /2 for 
2 < j < N. Like in the previous schemes, RUDD and 
RUDDp are differed by the constant phase and the phase 
alternation in 6. 



CPMG and CPMGp schemes involve periodic spin fiips 
as shown in Figure [Ijc). The phase is set to x in 
CPMG, while it is alternated between x and —x in CP- 
MGp. CPMG and CPMGp have different performances 
depending on the initial states 19, 34, 35]. The total du- 
ration of the N-pulse CPMG is T = A^(2t + t^), where 2r 
is the delay between the tt pulses and r^r is the duration 
of each tt pulse. The same parameters N and T are used 
to compare CPMG with the following schemes. 



B. UDD and UDDp 



The pulse distributions for UDD and UDDp schemes 



are shown in Fig. 
but not periodic 
time instant 



[d) . Here the spin flips are symmetric 
J|. The j*'^ TT pulse is applied at the 
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where T is the total duration of the sequence and N is 
the total number of pulses. For a finite bandwidth case, 
with a TT pulse of duration t-^ , the delays tj are given by 
Ti = TN+1 = ii -T^/2, Tj ^ tj+i-tj~T^, for 2 < j < N. 
Like in the previous scheme, UDD and UDDp are differed 
by the constant phase and the phase alternation in (j). 



IV. EXPERIMENT 

The sample consists of crystallites of powdered hexam- 
ethylbenzene. At room temperature the entire molecule 
undergoes six fold hopping about the Cq axis of ben- 
zene ring Further the methyl group rapidly reorients 
about its C3 axis. Due to these motions, the intramolecu- 
lar dipolar interactions are averaged out. Intermolecular 
dipolar coupling is retained and each molecule acts as 
a point dipole [7]. Under free precession (no DD), this 
sample has a spin-spin relaxation constant of about 25 
fis and a spin-lattice relaxation constant of 1.7s. All the 
experiments are carried out on a Bruker 500 MHz spec- 
trometer at room temperature. 



A. DD on SQC 

First we describe the performance of various DD 
schemes on SQC. SQC was prepared by using an ini- 
tial {tt /2)y pulse on equilibrium longitudinal magnetiza- 
tion. As described in section and Fig. [TJ;, CPMG se- 
quences were constructed by periodic distribution of tt 
pulses in T — TT — T fashion. The minimum r in our ex- 
periments was set to 2 /j,s owing to the duty cycle limit 
of the probe coil. Also duration of tt pulse was found to 
be = 4.3 /iS at maximum allowed amplitude. Under 
these experimental conditions, the allowed values of N 
for UDD and RUDD are 1 to 7. For TV > 8, one obtains 
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FIG. 3: Performance of 7-pulse DD cycles on SQC. The 
graphs correspond to r = 2 /is (a), r = 4 /is (b), t = 8 fis (c), 
r = 16 /is (d), T = 32 us (e), and t — 64 fis (f). 



negative delays between the pulses. Therefore to study 
DD schemes for longer durations, we cycled these 7-pulse 
DD sequences. The 7-pulse CPMG has a total cycle time 
of T{t) = 7(2t + Ttt)- The results of these experiments 
are shown in Fig. |31 The graphs correspond to r values 
ranging from 2 /xs to 64 /zs. The corresponding T(t) val- 
ues are used to select the sampling points in no DD, as 
well as to construct other DD sequences. It is clear from 
these data that RUDD displays superior performance for 
shorter r values, while CPMG shows better performance 
for longer t values. We can also notice from these plots 
that the performance of RUDDp is better than CPMGp 
and UDDp which have almost same performance. UDD 
has better behavior than RUDDp, and for longer r values 
UDD and RUDD have same behavior. However, these 
performances in general may also be dependent on ini- 
tial states p^ . 



B. DD on MQC 

As described in section HIl and Fig. HJa-b), the scheme 
for studying DD on MQC involves preparation of MQC, 
evolution of MQC, followed by storing MQC onto lon- 
gitudinal spin-orders. A delay t/j = 5 ms allows spu- 
rious coherences to decay. The longitudinal spin order 
is then converted into SQC using a (7r/2)j, pulse, fol- 
lowed by a purge pulse px of duration 50 /zs. A 180 
degree phase alternation of the detection pulse, purge 
pulse, and the receiver is used to reduce artifacts aris- 
ing from receiver ringing pzj. For efhcient generation of 
MQC five cycles of 8-pulse sequence shown in Fig. [TJd 
was used in preparation and mixing periods, and the pa- 
rameter A was optimized to 2 /is. In our experiments 



FIG. 4: Multiple quantum spectra showing different coher- 
ence orders detected after inserting different delays (a) and 
RUDD sequences of same durations (b). 



the coherences of successive quantum numbers are sep- 
arated by Aw = 27r X 200 kHz. In order to separate 
a maximum of rimax coherences, the relative phase a 
between the preparation and mixing is incremented by 
Aa = 7r/nmax, where rimax was chosen to be 64. The 
corresponding increment in the evolution period is given 
by Ati = Aa/Aw. The signal intensities of the spectrum 
corresponding to these increments after cosine transform 
display strong peaks at even multiples of Aoj. Mean value 
of the signal intensities is made to zero to suppress strong 
zero-quantum peak. Fig. |4^ displays these even MQCs 
detected after inserting various delays, and Fig. HJ) dis- 
plays those detected after applying RUDDp sequences 
of same durations. The first spectrum corresponding to 
no-delay is same in both cases, in which one can easily 
observe MQCs of order up to 14. Other spectra in (b) 
were obtained by RUDDp sequences constructed with in- 
creasing number of pulses, i.e.,7V=l,2,-- - ,7. Under no 
DD (Fig. |4^), the intensities decay monotonically with 
delays, while under RUDDp (Fig. |3)d) the dependence of 
intensities is oscillatory w.r.t. N. Similar behavior was 
earlier observed in a two-qubit liquid state NMR system 
[3l| . The spectra in (b) at odd N clearly show better 
intensities compared to the corresponding spectra in (a) . 
Comparisions of performance of different DD schemes for 
preserving MQCs of various orders are described in the 
following. 

The intensities of MQCs of even orders between 2 and 
8 w.r.t. size of various DD schemes are plotted in Fig. 
[5l Although we can observe MQCs of higher orders than 
8 for short durations, for longer durations their intensity 
measurements are not reliable due to noise. The first 
data point in each data set corresponds to no DD, and the 
rest correspond to different size of the DD sequence with 

= 1, 2, • • • , 7. As observed in Fig. Hb, we see the oscil- 
latory behavior of each MQC under various DD schemes. 
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FIG. 5: Performance of various DD schemes in preserving 
MQCs of order 2 (a), 4 (b), 6 (c), and 8 (d). Each data set 
has 8 points, in which the first point corresponds to no DD, 
and the rest correspond to different size of the DD sequences 
with = 1,2, • • • ,7. 
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FIG. 6: Performance of various DD schemes in preserving 
MQCs of order 2 (a), 4 (b), 6 (c), and 8 (d). The intensities 
were obtained from spectra detected after applying up to a 
majcimum of 3 cycles of 7-pulse DD schemes. 



The intensities of MQCs of even orders between 2 and 
8 for different cycles of 7-pulse DD schemes are plotted in 
Fig. ini The first data point in each case corresponds to 
no DD. The fast decay of magnetization under no DD al- 
lowed to detect intensities corresponding to a duration of 
only one cycle, while for RUDD and RUDDp, intensities 
up to 3 cycles could be detected. 

V. CONCLUSIONS 

We studied the performance of various DD schemes on 
nuclear spins with long-range interactions using a solid 
state NMR system. First we applied these DD schemes 
on a single quantum coherence. The experiments were 
carried out for various 7-pulse DD schemes and for dif- 
ferent delays (r) between the tt pulses. The results clearly 
show that all the DD schemes are able to preserve the sin- 
gle quantum coherence for longer durations of time com- 
pared to no DD. However, for small delays between the tt 
pulses, RUDD showed the best performance. For longer 
delays between the tt pulses, CPMG was better. Then 
we prepared MQCs of even orders using multiple cycles 
of the well known 8-pulse sequence implementing a two- 
quantum average Hamiltonian. The MQCs so prepared 
could be detected using standard spin-counting type ex- 
periments. We found the Various DD schemes were in- 
serted after the preparation of MQCs. We studied the 
performance of these DD sequences with different sizes. 
The intensity behavior under all the DD sequences were 
oscillatory. DD sequences with odd number of tt pulses 
showed much better performance than those with even 
number of tt pulses. RUDDp sequence showed the best 
performance over all other sequences. Since much of the 
theoretical and experimental work on dynamical decou- 
pling is limited to one or two qubits, we hope that the 
present work will shed light on the behavior of dynamical 
decoupling in systems with a large number of correlated 
qubits. It may also be interesting to study the behav- 
ior of odd quantum coherences under DD. It is possible 
to excite MQCs of all orders (both even and odd) si- 
multaneously using pulse sequences implementing single 
quantum average Hamiltonians [36j . 



But all the DD schemes display an overall improvement 
w.r.t. no DD. However it can be noticed that RUDDp has 
significantly better performance than all other schemes, 
even for higher order coherences. Surprisingly, unlike the 
single-quantum case, wherein RUDD displayed the best 
performance, in multiple-quantum case RUDDp is the 
best scheme. 
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